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Based on a general approach to binary systems we show that in low density region 
asymmetric nuclear matter (ANM) is unstable only against isoscalarlike fluctu¬ 
ations. The physical meaning of the thermodynamical chemical and mechanical 
instabilities is related to the inequality relations verified by the strength of inter¬ 
action among different components. Relevance of these results in bulk and neck 
fragmentation is discussed. 


1. Instabilities and fluctuations in ANM 

A binary system, including ANM, manifest a richer thermodynamical be¬ 
haviour as a consequence of the new degree of freedom. The phase tran¬ 
sitions are more complex because they have to accommodate one more 
conservation law. Also a variety of thermal fluctuations with different com¬ 
ponents composition can develop in the system. 

In particular, in symmetric nuclear matter (SNM) one encounters two 
kinds of density fluctuations: (i) isoscalar, when the densities of the two 
components oscillate in phase with equal amplitude; (ii) isovector when the 
two densities fluctuate still with equal amplitude but out of phase. These 
are the normal modes of the system in the sense that the variation of free 
energy density up to second order reduces to a normal quadratic form when 
is expressed in terms of them. Mechanical instability is associated with 
instability against isoscalar fluctuations leading to cluster formation while 
chemical instability is related to instability against isovector fluctuations, 
leading to species separation. 
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A general analysis, for ANM, can be performed in the framework of 
the Fermi liquid theory We limit ourselves to monopolar deformation 
for proton and neutron Fermi seas iz®, iz® and consider here momentum 
independent interactions such that {qi = n,p) is the only non-zero 

Landau parameter Then the thermodynamical stability condition at T = 0 
is obtained as: 

SH - UpSpp - PnSpn = > 0 (1) 

where H is the energy density and 

a = iVp(0)(l + FP^’) ; 6 = iV„(0)(l + Fo"") ; 

c=fV„(0)FP” + 7Vp(0)F7. (2) 


with Nq the single-particle level density at the Fermi energy We diagonalize 
the form of Eq.(l) by introducing the following transformation: 

u = cos/3 Up + sinf} u'^, 

V = -sinfi Up + cos/S u^, (3) 


where a mixing angle 0 </ 3 < 7 r/ 2 is defined by 

c _ /V„( 0 )FP" + /Vp( 0 )F 7 


tg 2/3 = 


a — b 


Np{0){l + PPP) - /V„(0)(1 + F™) ■ 


Then Eq.(l) takes the form 

SH — PpSpp — PnSpn = Xu^ + Yv^ 


( 4 ) 

( 5 ) 


where 

X = a -b 6 + sign{c)^/{a - by + + ^»( 0 ) ( i ps^ ) ( 5 ) 

and similar 

Y=\{a + b- sign{c)^{a-bf+c^ ) ^ + 1 + p^^ ). ( 7 ) 

Thus with the rotation (3) we separate the total variation Eq.(l) into 
two independent contributions, the new ” normal” modes, characterized by 
a ” mixing angle” /3, which depends on the density of states and the interac¬ 
tion. since cos/3, szn/3 are both positive, we interpret u- and n-variations as 
new independent zsosca/ar-like and isovector-\\\e fluctuations appropriate 
for asymmetric systems. The proton and neutron densities will fluctuate in 
phase for isoscalar-like variations and out of phase for isovector-like vari¬ 
ations, see Eq.(3) and they reduces to the usual isoscalar and isovector 
modes for SNM. 
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From Eq.(5) we see that thermodynamical stability requires X > 0 and 
Y > 0. Equivalently, the following conditions have to be fulfilled: 


1 + Fgg > 0 and 1 + F^g > 0, 


( 8 ) 


They represent Pomeranchuk stability conditions extended to asymmetric 
binary systems. 

The new stability conditions, Eq.(8), are equivalent to mechanical and 
chemical stability of a thermodynamical state, i.e. 




> 0 and 


T,y 


f 

\ J T,p 


> 0 


(9) 


where P is the pressure and y the proton fraction, see ® . Such general 
analysis leads to the conclusion that in the low energy region, where c < 0 
as predicted by all effective interactions, the asymmetric nuclear matter 
is unstable only against isoscalarlike modes. This is shown for a Skyrme- 
like interaction ® in Eigure 1, but this effect is expected to by very robust, 
present for all interactions. 



0.00 0.05 0.10 0.15 

Figure 1. Density dependence of general¬ 
ized Landau parameters for y=0.1 



Figure 2. The unstable normal modes for 
liquid phase in the spinodal region. 


The direction of unstable mode in the plane — pp in the liquid phase 
is indicated by the arrows in the Eigure 2. The arrow picks toward less 
asymmetric liquid phase: the isospin distillation take place. In this figure 
the length of the arrows was drawn to correspond to the same proton density 
perturbation. Therefore it indicate the ’’response” of neutrons and their 
relative sizes allows a comparison of the effect throughout spinodal region. 
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In the gas phase the asymmetry is larger then the initial value, the arrows 
being oriented in opposite directions. The tilting angle is a measure of 
isospin distillation which will depend on the slope of the symmetry energy 
in the low energy region. Therefore the isospin content of the fragments 
will be an interesting observable to constrain the isovector part of the EOS, 
see the next sections. 

The physical meaning of the thermodynamical chemical and mechani¬ 
cal instabilities is connected to the inequality relations existing between the 
interactions among different species ®. Indeed the system will show up as a 
chemical instability if {—ta — b/t) < c < —2\/^ and as a mechanical insta¬ 
bility if c < {—ta — b/t) < —2Vab, where t = jv**(oj • regions were 
each of these conditions are verified are indicated in Figure 3 as chemical 
and mechanical instability respectively. 



Q jji I I I I I I I I I.till i t.I. Pi I I I I I il ..I I I I I I I il 

0.00 0.02 0.04 0,06 0.08 0.10 0,00 0.02 0.04 0.06 0.08 0.10 0,00 0.02 0.04 0.06 0.08 0,10 

p(fm”^) P(fm“^) 


Figure 3. Spinodal boundary of asymmetric nuclear matter(open circles) and mechan¬ 
ical instability boundary (crosses) for three proton fractions:(a) y = 0.5,(a) y = 0.25,(a) 

?; = o.i 


2. Kinetic of Phase Transition in ANM 

Contrary to the situation regarding the phase transition in macroscopic 
systems were the observation time scale is much larger than the time scale 
of microscopic process that leads to drops (bubbles) formation (except for 
the critical point where are required days because of slowing dawn phe¬ 
nomenon) for the heavy ions collision the reaction time can be comparable 
to the fragment formation time. The violent collision may quench the sys¬ 
tem in the instability region of phase diagram. The kinetic mechanism 
responsible for the phase transition is the spinodal decomposition. The un¬ 
stable fluctuations are amplified and growth exponentially in the first stages 
leading to higher density domains. In asymmetric nuclear matter they are 
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especially interesting because of isoscalar like character of the unstable mode 
and presence of isospin distillation. 

We have studied the spinodal decomposition for ANM in a box with side 
L = 24: fm for several values of the initial asymmetry and initial density at 
temperature T = 5 MeV. Several kinds of initial perturbation were created 
automatically due to the random choice of particle positions. Results for 
the initial asymmetry / = 0.5 and initial conditions corresponding to chem¬ 
ical instability region {pmit = 0.09 fm~^) are reported in Figure 4 (a) and 
(b) respectively. It is shown the time evolution of neutron (thick histogram) 
and proton (thin histogram) abundance Figure 4 (a) and asymmetry ,Fig¬ 
ure 4 (b), in various density bins. The dashed line respectively shows the 
initial uniform density value p ~ O.Gpo and the initial asymmetry I = 0.5. 
The in-phase driving to higher density for neutrons and protons is evident 
indicating the isoscalarlike character of the instability. The consequence is 
the fragment formation. However the clustering mechanism is accompanied 
by isospin distillation leading to very different asymmetries, lower (higher) 
in the liquid (gas) phases. 

We want to stress that an identical qualitative behaviour is observed 
if we start from the mechanical instability region in agreement to the 
discussion of previous section. 


t=0 fm/c 


t=100fm/c t=200 fm/c t=250 fm/c t=300 fm/c 
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Figure 4. Time dependence of the proton and neutron numbers (a) and asymmetry (b) 
in different density bins. 
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3. Isospin dynamics from bulk to neck fragmentation 

We follow the reaction dynamics for the collision of the system 

at 50AMeV for two impact parameters b = 2/to and b = 6 /to respectively. 
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We consider an asystifF EOS corresponding to a linear dependence with the 
density of the mean-field symmetry term. 

For a semicentral collision, b = 2/m, the reaction mechanism corre¬ 
sponds to bulk fragmentation. After a hrst compression phase (until about 
40-50fm/c) a fast expansion phase follows (until 110-120fm/c). Then during 
the fragmentation stage the system will break up and the fragment forma¬ 
tion process take place up to the freeze-out time (around 260-280fm/c). 
These three main stages are characterized by specific features of the isospin 
dynamics since the system explores different density regions. 
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Figure 5. Time dependence of the asymmetry in gas and liquid phases for b=2fm (a) 
and b=6fm (b) as well as the asymmetry of the primary fragments at freeze-out time for 
b=2fm (c) and b=6fm (d). 


In particular we focused on the time evolution of the asymmetry of 
the ’’central” gas (solid line and squares), and total gas (dashed-bsquares), 
as well as of the ’’centraf’ liquid (solid-bcircles) and IMF (clusters with 
3 < Z < 23 , stars), 5 (a). By ’’central” we mean a region having a linear 
dimension of the order of 20/m corresponding to the active volume in which 
the fragmentation take place. 

The spinodal mechanism is well evidenced by the inversion in the trend 
for both liquid and gas phase corresponding to isospin distillation. The 
isospin content of primary fragments is below the initial value but also 
clearly lower than the value at the beginning of the phase transition, due 
to isospin distillation, see Figure 5 (a) and (c). 

For b=6fm we observe a quite different behaviour. Now in the overlap 
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region a neck structure is developing. During the interaction time (from 
about 80-120fm/c) it heats and expands but remains in contact with the 
denser and colder region of projectile-like (PLF) and target-like (TLF). 
Now the surface instabilities of a cylindrical shape neck region and the fast 
leading motion of the PLF and TLF will play an important role in the 
mechanism of IMF production. A correspondingly different dynamics of 
the isospin is evident from the Figure 5 (b). A larger asymmetry of IMF, 
Figure 5 (d), suggest a neck region richer (poorer) in neutrons (protons). 
This can be related to proton and neutron migrations between more diluted 
neck region and denser PLF and TLF as dictated by the chemical potential 
gradients of the two species. 


4. Conclusion 

In the previous section in an approach based on a Fermi liquid theory 
for two components systems we analyze several properties of ANM in low 
density region of phase diagram. It is concluded that: 

-in the low density region, of interest for nuclear liquid-gas phase tran¬ 
sition, the system can be characterized by a unique spinodal region, see 
also defined by the instability against isoscalarlike fluctuations; inside 
this we can identify a mechanical instability region and a chemical one; 

-the physical meaning of chemical and mechanical instabilities is con¬ 
nected to the inequality relations which are satisfied by interactions among 
different components; when this inequality changes the sign the chemical 
instability transform in mechanical one; 

-during the time development of the spinodal decomposition in ANM 
the fragment formation accompanied by the isospin distillation is observed. 
The qualitative behaviour is unchanged when we move the initial conditions 
from chemical instability region to mechanical instability one; 

We have also made a connection of these features with isospin transport 
properties in microscopic simulations of fragmentation reactions of charge 
asymmetric ions in the medium energy range. We observed a different 
isospin dynamics when we pass from bulk fragmentation, for semi-central 
collisions to neck fragmentation, in semi-peripheral collision. 

We also mention that a dependence on the symmetry term of EOS was 
evidenced for the observable showed above as well as for some other as is 
discussed in more detail in the refs. 
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